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ABSTRACT 

We present Ha and [O III] (5007 A) images of the nebula powered by the jet of the black 
hole candidate and microquasar Cygnus X— 1, observed with the 2.5m Isaac Newton 
Telescope (INT). The ring-like structure is luminous in [O III] and there exists a thin 
outer shell with a high [O III] / Ha flux ratio. This outer shell probably originates in the 
collisionally excited atoms close to the front of the bow shock. Its presence indicates 
that the gas is shock excited as opposed to photoionised, supporting the jet-powered 
scenario. The shock velocity was previously constrained at 20 < v s < 360 km s ; here 
we show that v s > 100 km s _1 (la confidence) based on a comparison of the observed 
[O III] / Ha ratio in the bow shock with a number of radiative shock models. From 
this we further constrain the time-averaged power of the jet: Pj c t = (4-14) x 10 36 



s . The Ha flux behind the shock front is typically 4 x 10 erg s 



erg 

1 cm~ 2 arcsec~ 2 , 

and we estimate an upper limit of ~ 8 x 10 -15 erg s _1 cm~ 2 arcsec -2 (3<r) to the 
optical (i?-band) continuum flux of the nebula. The inferred age of the structure is 
similar to the time Cyg X-l has been close to a bright H II region (due to the proper 
motion of the binary), indicating a dense local medium is required to form the shock 
wave. In addition, we search a > 1 degree 2 field of view to the south of Cyg X-l in 
Ha (provided by the INT Photometric Ha Survey of the Northern Galactic Plane; 
IPHAS) for evidence of the counter jet interacting with the surrounding medium. 
Two candidate regions are identified, whose possible association with the jet could be 
confirmed with follow-up observations in [S II] and deeper observations in [O III] and 
radio. 

Key words: accretion, accretion discs, black hole physics, ISM: jets and outflows, 
shock waves, X-rays: binaries, stars: individual: Cygnus X-l 



1 INTRODUCTION 

Galactic black holes release an unknown fraction of their 
accreting matter and energy in the form of collimated out- 
flows (or jet s) that travel into the surrounding medium (for 
a review see lFend"eii 2006). These accreting black hole can- 
didate X-ray binaries (BHXBs) spend a large fraction of 
their lifetimes in a 'hard' X-ray st ate (|Homan fc Bellonil 
120051 : iMcClintock fc Remillardll2006l ) which is ubiquitously 
associated with the formation of pa rtially self-absorbed syn- 
chrot ron emitting radio jets (e.g. iFender, Belloni fc Gallol 
120041 ). While the presence of these jets is commonly in- 
ferred from a flat radio-through-infrared spectral com- 
pon ent, radio observations of the nearby (d = 2.1±0.1 
kpci lMassev. Johnson fc Degioia-Eastwo od 1995) high-mass 
BHXB Cygnus X-l have directly resolved a collimated jet 
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~30 AU in length (|Stirling et alj|200ll ). whilst in this X-ray 
state. A transient radio jet has also been observed of length 
~140 AU th at was launched du ring a period of X-ray state 
transitions (|Fender et al.l |2006| ). Transient jets with high 
Lorentz factors are sometimes seen w hen a BHXB leaves 
the hard state (e.g. IFender et al.ll20o3 ). 

The power carried by the self-absorbed, steady hard 
state jets when averaged over their lifetime stands as one of 
the key questions in our understanding of BHXB accretion, 
and possibly for black hole accretion on all mass scales. In 
recent years it has become apparent that a significant frac- 
tion of the liberated acc retion power is ejected as outflo ws in 
hard state BHXBs (e.g. lGallo. Fender fc Poolevll2003l ). and 
it now seems likely tha t in hard X-ray states the acc retion 
flow is 'jet-dominated' l|Fender. Gallo fc Jonker| [2003). with 
more power in the jets than in X-rays. Clearly, measuring 
as accurately as possible the jet power is key to understand- 
ing both the overall physics of the accretion process and the 
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Figure 1. Stacked WFC images of the field of Cyg X-l. Top: 2400 
sec in Ha (6568A); Centre: 2400 sec in [O III] (5007A); Lower: 
1800 sec [O III] (5007A) of the field south of Cyg X-l. Cyg X-l 
is indicated by an arrow in each panel. In the top and middle 
panels, the H II region Sh2— 101 is prominent on the eastern CCD 
chips and the Cyg X— 1 shell nebula is visible to the north west of 
the system. 



matter and energy input from BHXBs into the interstellar 
medium (ISM). 

Attempts at estimating the jet power content from 
core radio luminosities are riddled with assumptions about 
its spectrum and radiativ e efficiency, the latter of which 
is poorly constrained (e.g . Oglev et all |2000| ; iFenderJ l200ll ; 
iHoman et all 120051 ; iHeind l2006h The jet power estimated 
in this way is highly sensitive to the location of the high- 
frequency break of the flat (spectral index a « where 
F v oc v a ) optically thick part of the jet spectrum, as the 
radiative power is dominat ed by the higher energy photons 
l|Blandford fc Konigl|[l979h . This break is likely to be at ~ 
2/jm, at least for BHXBs at high luminosities (Lx 10 36 erg 
i n the hard state (ICorbel fc Fenderll2002l ; iNowak et ail 



l2005l ; lRussell et al.|[2006f K 

The jet power may alternatively be constrained by 
analysing its interaction with the surrounding medium, 
without requiring prior knowledge of the jet spectrum and 
radiative efficiency. Radio lobes associated with jets from 
Active Galactic Nuclei (AGN) are commonly used as accu- 
rate calorimeter s of the power x lifetime product of the jets 
(Burbidgc 1959), a method only very recently applied to jets 
from stellar mass black holes. Radio lobes have been identi- 
fied and associated with an increasing number of BHXBs 
dMirabel et all 1 19921 ; iRodrfguez et al.1 Il992l ; ICorbel et al.l 
l2002h and even a couple of neutron star X-ray bina- 
ries dFomalont. Geldzahler fc Brad shaw 2001; T udose et al.l 
l200fl) . 

Recently, deep radio observations of the field of Cyg 
X-l resulted in the discovery of a shell-like structure 
which is aligned w ith the aforementioned resolved radio jet 
llGallo et aU 20051 '). This jet-blown nebula with a diameter 
of ~ 5 pc is interpreted as the result of a strong shock that 
develops at the location where the collimated jet impacts 
on the ambient ISM. The nebula was subsequently observed 
at optical wavelengths by the Isaac Newton Telescope (INT) 
Wide Field Camera (WFC) and is clearly visible in Ha (6568 
A) line emission, and coincident with the radio shell. The 
Ha flux density F„ > 0.02 mjy arcsec -2 is > 20 times the 
measured radio flux density, inconsistent with optically thin 
synchrotron emission between the t wo spectral regime s and 
instead favouring a thermal plasma (|Gallo et al.ll2005l ). 

Models of X-ray binary jet-ISM interactions predict a 
shell of shocked compressed ISM emitting bremsstrahlung 
radiation containi ng a bubble of rel ativistic synchrotron- 
emitting plasma l|Kaiser et al.l 120041 ). As with any inter- 
stellar shock wave where the preshock gas is at least par- 
tially ionised by the approaching radiation field, line emis- 
sion f rom recom bination of the shocked gas is also expected 
(e.g. ICoxlll972h . The shell, which is essentially a radiative 
shock wave, consists of a bow shock front where the gas is 
in transition to a higher temperature due to collision s and 
plasma instabilities (e.g. iMcKee fc Hollenbachl 11980!). The 
post-shocked gas then enters the optically thin 'cooling re- 
gion', where the overall ionisation level rises then falls as 
the g as radiatively cools to ~ 10 4 K (e.g. ICox fc Raymond! 
Il985h . Shock waves are commonly observed in supernova 
remnants (SNRs), and shocks created from outflows exist 
in Herbig-Haro objects, where emission lines are produced 
in the shock wave created by bipolar flow s from pre-main - 
sequence stars impacting the ISM (e.g. ISchwart3 Il983l ). 
Shock fronts associated with AGN jets interacting with the 
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Table 1. List of observations with the INT WFC. 



Run no. MJD Target RA & Dec. (centre of Airmass Filter Exposure 

CCD 4; J2000) time (sec) 



473434 


53653 


922 


Cyg X— 1 north 


19 


59 


38 


19 


+35 


18 


20 





1 


171 


[O III] 5007 A 


1200 


473435 


53653 


937 


Cyg X— 1 north 


19 


59 


39 


31 


+35 


18 


29 


6 


1 


231 


[O III] 5007 A 


1200 


473436 


53653 


956 


Cyg X-l north 


19 


59 


38 


50 


+35 


18 


20 





1 


310 


Harris V 


600 


473437 


53653 


963 


Cyg X-l north 


19 


59 


38 


49 


+35 


18 


20 





1 


380 


Ha 6568 A 


1200 


473438 


53653 


978 


Cyg X— 1 north 


19 


59 


39 


31 


+35 


18 


30 





1 


494 


Ha 6568 A 


1200 


473439 


53653 


996 


Cyg X— 1 north 


19 


59 


37 


68 


+35 


18 


10 





1 


620 


Harris R 


600 


473441 


53654 


013 


Cyg X-l south 


19 


59 


07 


99 


+34 


56 


59 


9 


1 


860 


[O III] 5007 A 


900 


473442 


53654 


025 


Cyg X-l south 


19 


59 


08 


80 


+34 


57 


09 


8 


2 


064 


[O III] 5007 A 


900 


473452 


53654 


106 


Landolt 94-171 


02 


53 


38 


80 


+00 


17 


17 


9 


1 


139 


Harris R 


10 


473454 


53654 


109 


Landolt 94-171 


02 


53 


38 


79 


+00 


17 


18 





1 


138 


Harris V 


10 


473457 


53654 


113 


Landolt 94-702 


02 


58 


13 


39 


+01 


10 


54 





1 


128 


Harris V 


10 


473458 


53654 


113 


Landolt 94-702 


02 


58 


13 


39 


+01 


10 


53 


9 


1 


128 


Harris R 


10 



The positions refer to the centre of CCD 4, which is the middle of the three eastern CCDs (Fig. 1). 10 bias frames and 5xHa, 7x [O III] 
5007 A, 6 X Harris _R-band and 5 X Harris V-band sky fiats were taken in the same observing night. 
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Figure 2. Image s of the Cyg X-l ring nebula in: (a) Ha; (b) [O III] ; (c) [O III] / Ha; (d) R-ban d - Ha and (e) radio 1.4 GHz 
l|Gallo et al.ll2005h . The Ha and [O III] images are de-reddened using Ay = 2.95 measured for Cyg X-l l|Wu et al.lll982li . flux calibrated 
and continuum-subtracted (i.e. Ha — R and [O III] —V). The [O III] / Ha flux ratio image is created from the non-continuum-subtracted 
[O III] and Ha frames. The Ha contribution has been subtracted from the i?-band image on the right, in which the nebula is not visible. 
Cyg X— 1 is the lower of the two bright stars near the bottom of each image and is marked by a cross in the radio image. 



intra-cluster gas (which is m uch hotter than the ISM) are 



also seen at X-rav energies (|Carilli, Perlev & Harris! 


1994; 


Kraft et all 120031; iForman et al.l 2005; Nulsen et all 


2005; 


Wilson, Smith & Youne 20061; Croston. Kraft & Hardcastle 



2006). In addition, infrared sources found symmet- 



ric about the position o f the X-ray binary GRS 
1915+105 dChatv et ai"]|200ll ) may be jet-ISM impact sites 
l|Kaiser et al.ll2004 ). 

By modelling the shell emission of the Cyg X- 
1 nebula as originating from radiatively shocked gas 
ilCastor. McCrav fc Weaver! 19751; iKaiser fc Alexanderl 19971 ; 
iHeinz. Reynolds fc Begelmanll 19981 ), the total power carried 
by the jet was estimated to be ~ 9 x 10 35 < Pj ct < 10 37 
erg s _1 which, when taking into account the power of the 
counter jet, is equal to ^ 0.06-1 times the bolometric X-ray 
luminosity of Cyg X-l (|Gallo et al.ll2005l ). These calcula- 
tions have le d to estimates of the matter content of the jet 
|Heindl2006T ); a similar technique to those applied to the 
jets of AGN, where their matter content are inferred from 
a combination of information from the core and lobe s (e.g. 
ICelotti fc Fabianll 19931 ; iDunn. Fabian fc Celottill2006l ). 



These jet power calculations are highly sensitive to the 
velocity of the gas in the shock front. From temperature con- 
straints and from the non-detection of an X-ray counterpart, 
this was estimated to be 20 < v s < 360 km s _1 . Emission 
line ratios of shock-heated gas can const rain its paramete rs, 
including the velocity of the shock (e.g. IOsterbrocklll9891 ). 

Here, we performed wide field imaging observations of 
the region of Cyg X-l in Ha and [O III] (5007 A) filters in 
order to (a) constrain the velocity and temperature of the 
shocked gas and hence the time-averaged power of the X- 
ray binary jet of Cyg X-l, (b) obtain for the first time flux- 
calibrated optical emission line measurements from a nebula 
powered by an X-ray binary jet, and (c) search for any jet- 
ISM interactions associated with the counter jet supposedly 
to the south of Cyg X-l. 



2 OBSERVATIONS AND DATA REDUCTION 

The field of Cyg X-l and of two Landolt standard stars were 
observed using the Wide Field Camera (WFC) on the 2.5m 
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Isaac Newton Telescope (INT) in narrowband and broad- 
band niters. The images were taken on 2005-10-10 as part of 
a Galactic survey of emission line nebulae powered by X-ray 
binary jets (Russell et al., in preparation). The conditions 
were cloudless with a seeing of ~1— 2 arcsec. The moon was 
set, making it possible to achieve a higher sig nal-to-noise 
ratio (S/N) than previous images presented in iGallo et al.l 
(2005). Table 1 lists the observations used for this work. 
The WFC consists of 4 CCDs of 2048x4100 imaging pixels, 
each of scale 0.333 arcsec pixel -1 . The field of view is 34x34 
arcmin. 

Data reduction wa s performed using the pipeline pack- 
age THELI (details in lErben et all I2005T ). After manually 
separating the data into type (bias/flat/science) and filter, 
THELI de-biased and flat-fielded the science frames. The 
package read the exposure time and RA & Dec from the 
FITS file headers (it is currently configured for reducing 
data from ~ 30 instruments including the INT WFC) and 
matched several hundred stars in each exposure with those 
in the online MAST Guide Star Catalog (release 2.2). THELI 
then aligned and stacked (for each position and filter) the 
images, and position-calibrated mosaics were created, with 
flux in counts s - pixel - . The resulting Ha and [O III] mo- 
saics of the Cyg X-l field are displayed in Fig. 1. 



3 FLUX CALIBRATION 

The Landolt standard stars 94-171 (V = 12.659; R = 
12.179) and 94-702 (V = 11.594; R = 10.838) were observed 
at low airmass (see Table 1). We performed aperture pho- 
tometry of the standards in the V and 7?-band reduced im- 
ages with PHOTOM and accounted for airmass-d ependent 
atmospheric extinction according to iKind l|l985h . The re- 
sulting conversion between intrinsic flux density (at airmass 
= 0) and counts s -1 differed between the two standards by 
a factor 0.5 and 6.5 percent for V and R, respectively. The 
Ha filter is located within the passband of the Harris R fil- 
ter, and [O III] is located within the Harris V filter. For a 
flat spectrum source, the WFC i?-band filter collects 15.742 
times more flux than the Ha filter, and likewise V collects 
11.127 times more flux than [O III]. By multiplying by these 
factors we obtain 7.140 x 10 -3 mjy (intrinsic) per counts 
s -1 for Ha and 5.637 x 10 -3 mjy (intrinsic) per counts s -1 
for [O III]. 

As a check of our calibration, we measured the fluxes 
of 10 stars in our reduced Cyg X-l Ha image, that are 
listed in the INT Photo metric Ha Surve y of the Northern 
Galactic Plane (IPHAS; lDrew et alj|2005h catalogue (which 
uses the same telescope, instrument and filter). The ratio 
of the Ha flux of each star measured here to that listed in 
the catalogue is 1.06±0.39 (lc). We can therefore assume 
our flux calibration is accurate, and we adopt a la error for 
each flux measurement F (Ha and [O III]), of 0.14 dex (i.e. 
log F±0.14). 

We created flux-calibrated continuum-subtracted Ha 
and [O III] images of the region of Cyg X-l (Fig. 2a, 2b) 
using IMUTIL and IMMATCH in IRAF . The resulting im- 
ages are corrected for interstella r extinction tow ards Cyg 
X-l adopting Ay = 2.95 ± 0.21 l|Wu et alj|l983 ). The ex- 
tinction towards Cyg X-l and its nebula may vary over 
the field of view. The total neutral hydrogen column den- 
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Figure 3. Total Galactic neutral hydrogen column density 
JVh in units of 10 21 cm -2 at positions close to Cyg X— 1 
iHartmann fc Burton1ll997l) . The box indicates the field of view 
(FOV) containing Cyg X-l and its nebula. 



sity Ah through the Galaxy at various positions in a ~ 4 
deg 2 area around Cyg X-l l|Hartmann fe Burtonlll997r i is 
indicated in Fig. 3. Ah appears to vary smoothly over 
the field, and the four closest points to Cyg X-l yield 
A H = (6.55T0.94) x 10 21 cm -2 in the region containing Cyg 
X-l, or A v = 3. 66±0. 53 a dopting A H = 1.79 X 10 21 cm -2 Ay 
l|Predehl fc SchmittfL 995). If the extinction Ay towards Cyg 
X-l (i.e. not through the whole Galaxy in this direction) 
also varies by ±0.53 over the much smaller field of view of 
the nebula (a conservative case), an error is introduced of 
0.18 dex and 0.25 dex to the Ha and [O III] flux measure- 
ments, respectively (no distance estimate is required). The 
measurements of [O III] / Ha however only suffer a 0.06 dex 
uncertainty due to the different extinctions suffered at the 
two wavelengths ([O III] at 5007 A and Ha at 6568A; this 
concerns only the contribution to the error budget of the 
extinction) . 

Propagating the above two sources of error, we arrive 

at: 

A(logi F Sa ) = 0.23 log 10 F Sa (1) 
A(log 10 F [o ni]) = 0.29 log 10 F [o m] (2) 
A{log 10 {F [o nq/Fna)) = 0.15 log 10 (F [o ui]/F Ha ) (3) 



[O III] / Ha and i?-band — Ha images were also cre- 
ated (Fig. 2c, 2d). The [O III] / Ha image represents the 
ratio of the [O III] and Ha fluxes (in de-reddened erg 
s -1 cm -2 arcsec -2 ) and the Ha contribution is subtracted 
from the continuum in the J?-band — Ha image. The 1.4 
GHz radio image of IGallo et ail l|2005l ) in shown in Fig. 2e 
for comparison. 



4 RESULTS 

The Ha and [O III] fluxes measured in a number of positions 
in the shell of the nebula are plotted in Fig. 4. It seems 
that the [O III] / Ha ratio increases with position going 
clockwise around the shell, from values of F\o nil / Fu a ~ 
0.2 on the eastern side to ~ 1.2 on the western side of the 
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Figure 4. Intrinsic de-reddened continuum-subtracted Ha and 
[O III] fluxes at various positions around the shell, versus the 
position angle from Cyg X—l. See equations 1-3 for the errors as- 
sociated with each flux measurement. At each position, the mean 
flux in a 164 arcsec 2 circular aperture containing no visible stars 
is plotted. The data on the left are from the eastern side of the 
nebula and the data on the right are from the western side. 



nebula. The fluxes of both emission lines are also higher on 
the western side. The line ratio is approximately constant 
in the filaments inside the nebula, with values of Fro m] / 
Fhc ~ 0.7. The H II region Sh2-101 and the 'stream' of 
diffuse emission joining the region to the Cyg X-l nebula, 
possess ratios of order F[o n i] / Fn a ~ 0-0.2, with one area 
(RA 19 59 45, Dec +35 18 48) having a higher value, F [Q IU] 
I Fna ~ 0.6. This area resembles a shell-like feature close 
to a bright star HD 227018 at the centre of the H II region, 
and could be a bow shock associated with the motion of 
this star (which is 12 milli-arcsec per year in the direction o f 
the observed shell in the H II region; IPerrvman et al.lll997t ). 
We measure the flux in the H II region to be Fu a ~ 5- 



15 x 10" 14 erg s 



_1 cm -2 arcsec -2 ; F [Q nl] ~ 0.5-5 x 10 -14 
erg s -1 cm -2 arcsec -2 . 

The large quantity of stars in the line of sight of the 
shell prevent us from accurately calculating the total Ha 
and [O III] luminosity of the structure. However, we can 
conservatively estimate the area of the emitting gas to be 
between 3 arcmin 2 (approximating the nebula to a 15 x 700 
arcsec 2 area) and 46 arcmin 2 (the area of a circle with radius 
4 arcmin) and take the mean flux of the nebula (Fig. 4) to 



constrain its luminosity. At a distance of 2.1 kpc, the Ha and 
[O III] luminosities of the nebula are 1.8 x 10 34 < Lu a < 
2.8 x 10 35 erg s -1 and 1.3 x 10 34 < L [Q m] < 2.1 x 10 35 erg 
S - , respectively. 

In ten regions (of ~ 5 arcsec 2 areas) inside and out- 
side the Ha-emitting nebula, the i?-band de-reddened Ha- 
subtracted flux is consistent with zero, and the mean of 
the standard deviations in each region is 2.57 x 10 -15 erg 
s -1 cm -2 arcsec -2 . We therefore have a 3cr upper limit to the 
i?-band flux of 7.7 x 10 -15 erg s -1 cm -2 arcsec -2 , or a flux 
density of 0.0068 mjy arcsec -2 . Similarly, we obtain a 3<r up- 
per limit to the l^-band flux (from ten ~ 5 arcsec -2 regions) 
of 2.0 x 10 -14 erg s -1 cm -2 arcsec -2 , or 0.0196 mjy arcsec -2 . 
From the 7?-band upper limit, we conclude that the optical 
continuum flux de nsity is < 7 times the radio flux density 
l|Gallo et al.ll2005h . yielding a radio-optical spectral index 
a < 0.15. 

The broadband spectrum is consistent with 
bremsstrahlung emission with a flat radio-optical con- 
tinuum and optical emission lines which are at least one 
order of magnitude more luminous than the continuum. We 
cannot however definitely state that the radio emission is 
thermal in origin, as mixed thermal/non-thermal plasma 
is seen in some environments (e.g. the arcsec-scale jets 
of SS 433; iMigliari et "all 120021 ). A future measurement of 
the radio spectral index, which is difficult to achieve, is 
required to test whether the radio through optical spectrum 
is consistent with a thermal plasma or requires two (or 
more) components. Currently the morphology of the radio 
emission, being spatially coincident with the rather thin 
optical shell, does seem to favour the single thermal plasma 
model. 



4.1 Morphology and nature of the nebula 

The higher S/N in the Ha image compared to that in 
iGallo et"ail (|2005l ) is striking from the observed structure 
within the nebula (Fig. 2a). One would expect the most lu- 
minous areas of the nebula to occur in regions of the shock 
front where the line of sight is tangential to its surface and 
therefore should exist only at the apparent edge of the neb- 
ula if the structure is spherical. Filaments of gas are seen 
within the nebula and are unlikely to originate in unassoci- 
ated line-of-sight photoionised gas because of their filamen- 
tal morphology; the nebulosity is less diffuse than the gas 
in the photoionised H II region. The [O III] / Ha ratios of 
~ 0.7 in the filaments (see above) and morphology suggest 
either multiple shock fronts or an uneven non-spherical neb- 
ula. The latter is expected for a varying preshock gas density, 
which is clearly the case on inspection of the variation in Ha 
nebulosity close to the nebula. 

The [O III] / Ha image reveals two striking features: 
firstly the Ha nebulosity (which extends from the H II region 
Sh2-101 to the east of the nebula) possesses a low [O III] / 
Ha flux ratio compared to the Cyg X-l nebula, and secondly 
there exists a thin outer shell to the Cyg X-l nebula with a 
higher [O III] / Ha ratio. H II regions and photoionised gas 
in ge neral can possess a wide range of [O III] / Ha ratios 
(e.g. |johnson|[l953l: iBlair, Kirshner fc Cheva lier 1982]) but 
this morphology of an outer [O IIl]-emitting shell is only ex- 
pected from shock-excited gas. These two observations con- 
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Figure 5. A blow-up of part of Fig. 2c: [O III] / Ha image of 
the Cyg X— 1 ring nebula. We interpret the thin outer shell with a 
high [O III] / Ha flux ratio as possibly originating in the ionised 
atoms close to the front of the bow shock. North is up, east to 
the left. The emission line fluxes and [O III] / Ha ratio along the 
slice indicated in the image are shown in Fig. 6. 



firm that the gas in the Cyg X-l nebula must be ionised by 
a different source to the surrounding nebulosity. 

In Fig. 5 we show a blow-up of the north-west area of 
the nebula in [O III] / Ha; clearly showing this thin shell of 
enhanced [O III] / Ha ratio on the outer side of the ring. 
We have taken a 3.3 arcsec-wide slice through the image 
(indicated in Fig. 5) and plot the emission line flux and 
F[o in] / Fua flux ratio along this slice in Fig. 6. The upper 
panel of Fig. 6 shows that there is a sharp outer edge to 
the ring, and this edge is further out in [O III] than in Ha. 
This results in the thin outer shell of enhanced [O III] / Ha 
ratio that is visible in the lower panel of Fig. 6. This outer 
shell is ~ 6 arcsec thick and the rest of the ring is ~ 40 
arcsec thick, at least in the region of the slice. We note that 
the thin shell of enhanced [O III] / Ha ratio is not due to 
inaccurate alignment of the [O III] and Ha images as the 
stars are aligned and the thin shell exists on both sides of 
the nebula. 



4.2 Constraining the shock velocity and jet power 

We measured the F[o ni] / Fna flux ratio in seven 'slices' 
(typical thickness 12 pixels = 4 arcsec) of the [O IIl]-emitting 
shock front, and in four larger regions behind it. In the thin 
outer shell we obtain Fjo mi / Fn a = 1.36T0.28 and behind 
the shell, Fro in] / Fua = 0.89 ± 0.46. Taking into account 
the error from equation 3, F] hi] / Fn a = 1.3&ta % in the 
shock front and F] hi] / Fna = 0.891q 49 in the nebula 
behind the shock front. 

Both radiative and non-radiative shock waves are ex- 
pected to produce a mix of line (radiative recombina- 
tion) and continuum (d ue to bremsstrahlung emission) out- 
put (e.g. iKrolikl I1999T ). Non-radiative shock waves (de- 
fined as when the cooling time of the shocked gas ex- 
ceeds the age of the shock) produce Balmer optical emission 



a 
K 



q 

Oh 




20 40 60 80 100 
distance along slice from SE end (arcsec) 

Figure 6. The emission line flux (top panel) and the [O III] / 

Ha ratio (lower panel) along a slice through the nebula indicated 

in Fig. 5. The [O III] flux has been offset by —5 X 10 -15 erg 
— 1 —2 



2 for clarity. In the lower panel, region a = inside 
the ring nebula; c = the thin outer shell luminous 



'cm "arcsec 
the nebula; b = 

in [O III] / Ha; and d = outside the nebula. The mean error 
associated with each data set (equations 1—3) is indicated on the 
left in the panels. The spikes in the upper panel arc due to stars 
in the slice. 



lines only, because the hea vier atoms are not excited (e.g. 
iMcKee fc Hollenbachl [l980T ) . so the existence of the strong 
[O III] line confirms the radiative nature of the Cyg X-l neb- 
ula shock. In a radiative bow shock, the ambient ISM gas 
is first perturbed by the approaching radiation field. The 
ionisation of preshock hydrogen becomes progressively more 
significant in the velocity range 90 < v B < 120 km s~ x ; be- 
low 90 km s" 1 the preshock gas is essenti ally neutral and 
abov e 120 km s _1 it is entirely ionised fe.g. IShull fc McKeel 
1979). In addition, at velocities ~ 100 km s 1 the ionis- 
ing He II (304A) photons from the cooling region can dou- 
bly io nise the oxygen increasing t he intensity of the [O III] 
lines (|McKee fc Hollenbachl Il980fi . This leads to the high 



in all shock-wave models ( 


Raymond 19791; Shull & McKeel 


19791 Dooita et al. 19841; 


Binette. Dopita & Tuohv 1985; 


Cox & Raymond| |l985|; |Hartigan, Raymond & Hartmannl 


1987t iDopita fc Sutherland! 1 19951; 1996). The outer shell of 



the Cyg X-l nebula, bright in [O III] / Ha, may originate in 
the ionised atoms close to the front of the bow shock, and 
its presence confirms its association with the jet of Cyg X-l. 
Radiative shock models predict that relative optical 
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Cyg X-l (lo error range) 
Cox & Raymond 1985 model 
Binette et al. 1985 model 
Hartigan et al. 1987 model E 
Hartigan et al. 1987 model D 
Dopita Sutherland 1995 model 



100 200 300 400 

v s (km s' ) 



500 



600 



Figure 7. The [O III] / Ha ratio as a function of shock velocity 
predicted from five radiative shock models. The dashed rectan- 
gle represents the constrained range of v alues measured fo r the 
Cyg X—l shock front (velocity range from lGallo et alj|2005h . The 
INT [O III] filter includes both lines of the doublet; 4958. 9A and 
5006.9A, and the Ha filter includes the two [N II] lines at 6548.ll 
and 6583. 4A in addition to 6562. 8A Ha, within its bandpass. We 
therefore plot log(F [0 in] + [o 11T\/ F [N II]+Ha+[N II]) versus veloc- 
ity. 



line strengths of the shocked gas are highly sensitive to 
the temperature and velocity of the shock wave, and vary 
only weakly with the conditions i n the preshock gas suc h 
as the elemental abundances (e.g. ICox fe Raymond! 1 1985h . 
Temperatures and velocities of shocked matter, mainly in 
SNRs, are commonly inferred by comparing their observed 

optical or ultraviolet (UV) line ratios to these models (e.g. 

— ll 1 l ll l ll — — i 

Blair et al. 1991; Vancura et al. 1992; Lcvcnson et al. 1995; 



Bocchino et al.ll200d ; lMavromatakis et al.ll200lT ). The bright 



[O III] / Hcv shock front of the Cyg X-l nebula appears to 
be thicker on the western side of the nebula, and the Ha 
and [O III] surface brightness is also higher in this region 
(Fig. 4) . This implies the shock may be faster here, possibly 
due to a lower preshock density; it is the furthest area of the 
nebula from the H II region (or due to the proper motion of 
Cyg X-l; see below). 

The relation between the shock velocity and [O III] / 
Ha ratio in the shell is described for most of the published 



pre-ionisation ( 


Cox & Ravmondll 19851: iBinette et al.lll985l: 


Hartigan ct al. 


1987;lDoDita & Sutherland 1995). The reli- 



ability of these models can be assessed from their level of 
agreement, unless there are systematic uncertainties com- 
mon to all models (which seems not to be the case as the 
velocities predicted by the models are c onsistent with veloc- 
ities observed using other methods; e.g. iMichael et all 2000; 
iPun et all 120021, . In Fig. 7 we plot the aforementioned re- 
lation predicted from five of the most advanced radiative 
shock models. We overplot the empirical [O III] / Ha ratio 
found for the Cyg X-l shock front. Bearing in mind the pre- 
viously estimated range of shock velocities, 20 < v s < 360 
km s~\ we find that outside the range 100 < v B < 360 km 
s" 1 , no models successfully predict the observed [O III] / Ha 
measurement to within la. Subsequentl y, we infer a tem - 
perature of the gas in the shocked shell l|Gallo et al.ll2005h 
of 2.4 x 10 5 < T s < 3.1 x 10 6 K. The models in Fig. 7 



assume typical Galactic abundances and preshock particle 
densities within the range inferred for the Cyg X-l nebula; 
1 < no < 300 cm -3 (|Gallo et al.ll2005h . 

By constraining the velocity and temperature of the gas 
in the shock front, we are able to refine a number of parame- 
ters such_as_the particle density in the shell and power in the 
jet |Gallo et alj|2005h . Table 2 lists the parameters inferred 
for the nebula and jet. We adopt t he same methodolo gy to 
calculate the parameters as used in lGallo et al.l (|2005l ). The 
temperatu re-dependent ionisation fraction of the gas x is 
taken from Spitzcr (1978). At the temperatures inferred, the 
hydrogen in the shocked shell is fully ionised. We have im- 
proved the accuracy of the shock velocity, the time-averaged 
jet power and the jet lifetime from > 1 dex to ~ 0.5 dex. 
Assuming the power originates in the hard state jet which 
is s witched on for ~ 9 percent of the lifetime of Cyg X-l 
(e.g. iGallo et al.ll2005f) . the total power of both these jets is 
between 30 and 100 percent of the bolometr ic 0.1-200 keV 
X-ray luminosity (Lx ~ 3 x 10 37 erg s _1 ; |Pi Salvo et al.l 
l200lh . 

Assuming the shock velocity has been approximately 
constant over the lifetime of the jet, the jet lifetime t ~ 0.02- 
0.04 Myr is much shorter than the est imated age of the pro- 
genit or of the black hole, ~ 5-7 Myr llMirabel fc Rodriguesl 
2003). However, given that the proper motion of Cyg X- 
1 is ~ 8.7 mas y ear -1 in a direction ~ 130 degrees from 
the jet direction ^Stirling et all l200ll : iLestrade et~ai1ll999l : 
iMirabel fc Rodriguesll2003l l the system must have travelled 
between 2.5 and 9.1 arcmins during the lifetime of the jet. 
We speculate that the jet may have been present for the life- 
time of Cyg X-l, but only when the system moved into the 
ISM close to the H II region did the local density become 
large enough for a bow shock to form. In fact, the age of 
the nebula (if it is powered by the jet) cannot exceed ~ 0.04 
Myr because the position of Cyg X-l at the birth of the jet 
would be in front of the bow shock. This scenario may also 
explain the apparent higher [O III] / Ha ratio and velocity of 
the gas on the western side of the nebula compared with the 
eastern side: the western side is further from the jet origin 
as the Cyg X-l system moves. 



4.3 Ruling out alternative origins of the shell 

We have shown (Section 4.1) that the morphology of the Cyg 
X-l nebula indicates the gas is shock-excited and not pho- 
toionised; the structure is therefore not a planetary nebula 
or an H II region. The ne bula also cannot b e a SNR associ- 
ated with Cyg X-l itself (|Gallo et al.ll2005h . The proximity 
of the nebula to Cyg X-l, its alignment with the Cyg X-l 
jet and the agreement of its age inferred from its velocity 
and from the proper motion of Cyg X-l (Section 4.2) all 
point towards the jet-powered scenario. The one remaining 
alternative is a field SNR not associated with Cyg X-l; we 
her e explore this possib i lity. 

IXu. Zhang fc Hani l|2005al ) plot the diameter of 185 
Galactic SNRs against their 1 GHz surface brightness (Fig. 3 
of that paper). The radio surface brightness of the Cyg X-l 
nebula is Ei ghz « 3 x 10 ~ 23 W m~ 2 Hz" 1 sr" 1 (lower than 
any SNR in the sample of lXu et alJl2005aT ). and would have 
a diameter of > 100 pc if it were to lie close to the diameter- 
surface brightness relation for SNRs. The nebula has an an- 
gular diameter of ~ 7 arcmin, suggesting a distance to the 
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Table 2. Revised nebula and jet parameters (la confidence limits 
where quoted). 



Parameter 



Value 



Unchanged parameters (Gallo ct al. 2005): 

Monochromatic radio L v 1.4GHz ~ 10 18 erg s _1 Hz _1 
luminosity 

Cyg X— 1 - shell separation L fa 3.3 X 10 19 cm 

Jet inclination angle 9 fa 35° 

to line of sight 

Thickness of shell AR fa 1.6 X 10 18 cm 

Source unit volume V(L, AR) fa 4 X 10 53 cm 3 

Bremsstrahlung emissivity e„(V, L„ 1.4GHz) ~ 2.5 X 10 — 36 

for hydrogen gas erg s _1 Hz _1 cm -3 

Gaunt factor g fa 6 

Number density of 1 < no < 300 cm -3 

preshock gas 

Refined parameters: 



Ha flux of shell 


Fna ~ 


S.6 x 10 -15 




erg s 


— 1 cm — 2 arcsec — 2 


Ha luminosity of nebula 


1.8 X 10 


34 < L Ha < 




2.8 x 


10 35 erg s" 1 


[O III] flux of shell 


F[o in] 


fa 2.8 x 10 -15 




erg s 


— 1 cm — 2 arcsec -2 


[O III] luminosity of 


1.3 x 10 34 < L ro mi < 


nebula 


2.1 x 


10 35 erg s _1 


[O III] / Ha ratio 


F[o in] 


/ F Ha = 1.36±° ,% 


in shock front 






[O III] / Ha ratio 


F [0 III] 


1 F Ha = 0.89l°-j 4 


behind shock front 






Optical continuum flux 


Fv,OPT 


< 6.8ujy arcsec - 2 


density of shell 


(3a) 




Shock velocity 


100 < v 


3 < 360 km s" 1 


Shocked gas temperature 


2.4 x 10 


5 < T s < 3.1 x 10 6 K 


Number density of ionised 


56 < n e 


< 103 cm" 3 


particles in shell 






Ionisation fraction 


x fa 1 




Number density of total 


56 < rip 


< 103 cm -3 


particles in the shell 






Jet lifetime 


0.017 < 


t < 0.063 Myr 


Time-averaged jet power 


P Jct = (4-14) X 10 36 erg s" 1 


Hard state jet power 


Pjets — 


(9-30) X 10 36 erg s -1 


including both jets 






Outflow power / X-ray 


0.3 < /j 


ct/X < 1.0 



luminosity 



source of d>50 kpc if it is a SNR. However, the authors 
note that this surface brightness-diam eter relation could 
be sev erely biased by selection effects (|Xu. Zhang fc Hani 
2005b). Large, low surface-brightness SNRs may not have 
been detected in surveys due to con fusion with strong fore- 
ground and background sources (e.g. Green 2004). We there- 
fore cannot completely rule out a SNR origin to the nebula, 
but its alignment with the relativistic jet of Cyg X-l (for 
example) strongly suggests it is related to the system. 



4.4 The search for ISM interactions with the 
southern jet 

We searched a ~ 0.25 deg 2 field in [O III] south of Cyg X-l 
for a shell powered by the southern jet (lower panel of Fig. 
1). No shell was observed in the 1800 sec image, nor any 



visual nebulosity detected in the [O III] images that is con- 
sistent with the southern jet of Cyg X-l interacting with 
the ISM. The mass density of ISM may be much less to the 
south; the jet should travel essentially ballistically a nd un- 
detect ed until being decelerated by a denser region. iHeinj 
(2002) found that the rarity of jet-powered lobes in micro- 
quasars compared to AGN is likely due to the lower density 
environment. In addition, it is interesting to note that the X- 
ray-bright shocked e dges surrounding A GN radio lobes are 
sometimes one-sided |Kraft et al.ll2005T ). We searched a ~ 1 
deg 2 Ha mosaic from the IPHAS survey l|Drew et al.ll2005l . 
see Fig. 8). Two candidate 'hot spots' analogous to those 
associated with AGN jets were identified, and are expanded 
in the figure. 

Radio emission e xists in the region c lose to the northern 
hot spot according to lMartf et al.l (|l996t ). but it is uncertain 
whether it originates in point sources or extended emission. 
We inspecte d the southern reg ion of the 1.4 GHz radio image 
obtained bv lGallo et al.l ([2005) to search for diffuse emission. 
In a 1600x1600 arcsec region close to the more northern Ha 
hot spot (the southern hot spot is not within the radio field 
of view), the mean r.m.s. flux density is 4 x 10 -4 Jy beam -1 , 
yielding a 3<r upper limit to the (putative) radio hot spot in 
that region of 1.2 mjy beam -1 . Deeper radio observations 
are needed to test whether the candidate Ha hot spots have 
radio counterparts, as would be expected if they are jet-ISM 
interaction sites. In addition, observations in [S II] would test 
whether the emitting gas is shock-excited or photoionised, 
and deeper [O III] observations may constrain the velocity 
of any shocked gas. 



5 CONCLUSIONS 

We have imaged the shell-like nebula supposedly powered by 
the jet of Cyg X-l in narrowband optical emission line fil- 
ters. The nebula is bright in Ha and [O III] and the morphol- 
ogy and line ratios indicate that it is a radiative shock wave 
with a velocity v B > 100 km s - . We rule out alternative in- 
terpretations of the structure such as photoionised gas, but 
cannot completely rule out an unrelated SNR origin. How- 
ever, the morphology of the nebula and its alignment with 
the jet of Cyg X-l strongly suggest that it is indeed powered 
by this jet. We have probably isolated for the first time the 
thin, hot, compressed gas close to the bow shock front of a 
nebula powered by an X-ray binary jet. The time-averaged 
power of the jet must be P 3ct = (4-14) x 10 36 erg s -1 to 
power a shock wave of the inferred size, luminosity, velocity 
and temperature. If the steady hard state jet of Cyg X-l 
(that is active for ~ 90 percent of th e time) is providing this 
power (rather than the transient iet; lFender et al]|2006 ). the 
total energy of both jets in the hard state is 0.3-1.0 times 
the bolometric X-ray luminosity. This provides the strongest 
evidence so far that at an accretion rate of ~ 2 percent Ed- 
dington, Cyg X-l is in a state of approximate equiparti- 
tion between radiative and kinetic outp ut. Lower accretion 
rates should therefore be jet-dom inated l|Fender et alj|2003l ; 
iKording. Fender fc Migliaril [20 06) . The inferred age of the 
nebula is consistent with the time Cyg X-l has been close 
to the currently nearby H II region. 

No nebula associated with the southern jet of Cyg X-l 
is observed, although we detect two putative Ha 'hot spots' 
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left. RA and Dec are epoch J2000. The most extreme opening position angles of the resolved radio jet (both the steady low/hard state 
and transient jets) north of Cyg X-l; 17° < < 28° Stirling et al.ll200ll ; lFender et alj|2006t also projected to the south assuming the 
jets are anti-parallel) are shown. The north nebula and two areas of candidate ISM interactions with the southern jet of Cyg X— 1 are 
expanded. 



that require optical emission line and radi o follow-up to con - 
firm their nature. Optical and UV (e.g. iBlair et al.l Il99ll ) 
spectra of the northern shell will provide invaluable insights 
into the properties of the shocked gas, and hence the jet. 
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